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We present model results for the two-halo-neutron correlation functions, Cnn, for the dissociation 
process of light exotic nuclei modelled as two neutrons and a core. A minimum is predicted for 
Cnn as a function of the relative momentum of the two neutrons, pnn, due to the coherence of the 
neutrons in the halo and final state interaction. Studying the systems ^"^Be, ^^Li and ®He within this 
model, we show that the numerical asymptotic limit, Cnn —» 1, occurs only for pnn ^ 400 MeV/c, 
while such limit is reached for much lower values of pnn in an independent particle model as the 
one used in the analysis of recent experimental data. Our model is consistent with data once the 
experimental correlation function is appropriately normalized. 

PACS numbers: 25.10.-l-s, 21.45.+V, ll.80.Jy, 21.10.Gv 


The discovery of radioactive exotic weakly bound nu¬ 
clei, rich in neutrons or protons, beyond the drip line 
have brought a lot of interest in the nuclear structure 
and reactions of these unstable nuclei. The traditional 
nuclear models are unable to describe the long-range cor¬ 
relations between the nucleons of the halo. The large size 
of the halo has stimulated several experimental and the¬ 
oretical studies devoted to clarify the new aspects of the 
structure, stability and reaction of these nuclei, including 
astrophysical applications (see, for example, [ESI3)- 

The effect of the large spatial extension of the neutron 
halo, was probed in a recent fusion experiment of ®He 
with 238 jy target Q. It was observed a large reaction 
cross-section due to a direct 2n (n represents a neutron) 
transfer rather than an enhancement of the complete fu¬ 
sion cross-section. In view of the large size of the halo 
(r™'* = 5.9± 1.2 fm [^, r™® is the nn root-mean-square 
radius) which is comparable to the size of 238 itself, it 
is most likely that a correlated neutron-neutron (nn) is 
transferred to 238 jy -vvhile the "'^He is still far from the 
fusion barrier, favoring the 2n transfer process. 

A novel view of the complete fusion reaction for large 
halo nuclei also emerges it is most probable that the 
core approaches the fusion barrier of an excited target 
nucleus that has already absorbed the halo. This in¬ 
terpretation is in line with the concept of an absorptive 
many-body potential written in the relative coordinates 
of the core, two neutrons and the target [3|- 

It is evident the need of a deeper insight into the key 
aspects of the structure and dynamics of the halo in Bor- 
romean three-body systems (where all the subsystems 
are unbound), like ®He. In this respect, using intensity 
interferometry ^rfied with a new iterative technique, 
Marques et al. ll| recently probed the spatial config¬ 
uration of two-neutron halo systems and estimated the 
mean-square nn distances as well, considering the dissoci¬ 
ation of ®He, ^^Li and ^"^Be in the field of a heavy nucleus 
target. The spatial conheuration of ^^Li was also stud¬ 
ied by Petrascu et al.|^|l3- The nn correlation function, 


Cnm is extracted as a function of the relative momentum 
between the neutrons, Pnn- As the absolute normaliza¬ 
tion of the correlation function is an important piece of 
physical information, one should consider a consistent 
model for that. We note that, for the asymptotic nor¬ 
malization, in the ht to data, it has been considered a 
model where the halo neutrons are assumed to be inde¬ 
pendent El In our understanding, one should assume 
them as being emitted by a coherence source. 

In the present communication we report our results 
for the Cnn of Borromean three-body systems n — n — A, 
where A is the core mass number. We consider ®He, ^^Li 
and ^^Be within a description of two neutrons and a core 
forming Borromean systems. We show how to circum¬ 
vent the major difficulty of the relationship between the 
initial and final states due to the distorting effects of the 
reaction: we consider the nn final state interaction (FSI) 
and the three-body structure of the Borromean system. 
The nn FSI is shown to play a crucial role in distorting 
the relative motion of the neutrons and it is a source of in¬ 
terference effects which leads to an unexpected minimum 
of the Cnm as a function of Pnm pushing the asymp¬ 
totic behavior (where Cnn = 1) to larger momentum. 
By comparing our results with the available experimen¬ 
tal ones [1,11,13, we found that it is reasonable to expect 
a different normalization for the data, as will be shown. 

We use a three-body model in the limit of zero-range 
interaction, which retains the essential physics of the 
weakly bound and large two-neutron halo systems H, |1] . 
The interaction singularity is tamed in a renormalized 
zero-range model Il2l HI HI HI which is appropriate 
to study weakly bound three-body systems. The model 
is parameterized by minimal number of physical inputs, 
which are directly related to known observables: the two- 
neutron separation energy, S[2n) = —Ennc, the nn and 
neutron-core (n — A) s—wave scattering lengths (or the 
corresponding virtual or bound state energies). 

The neutrons of the halo have a large probability to be 
found outside the interaction range. Therefore the low- 
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energy properties of these halo neutrons are, to a large 
extend, model independent as long as few physical input 
scales are fixed. The model provides a good insight into 
the three-body structure of halo nuclei, even considering 
some of its limitations. It is restricted to s—wave two- 
body interactions, with small energies for the bound or 
virtual states. Even in this case, the three-body wave 
function for the valence neutrons presents configuration 
mixin g d ue to angular momentum recoupling [see also 
Refs. |l4 IT^l. We also note that, all the interaction 
effects, such as higher partial waves in the interaction 
and/or Pauli blocking effect are, to some extend, included 
in the model, as long as the experimental two- an three- 
body energies are supplied. In this perspective, it is of 
interest to extend the three-body continuum calculation 
of Cnm done in da for ®He, to analyze the actual data. 

The zero-range interaction itself does not have a physi¬ 
cal scale, then any calculated observable only depends on 
the two and three-body energy scales. The dependence 
on the three-body scale is a consequence of the Thomas 
collapse of the three-body system in the limit of a 
zero-range interaction which demands a three-body scale 
to supply the system with the physical information of the 
three-particle short distance configurations. Therefore, 
any dimensionless observable is a function of the ratios 
between the input energies, the mass ratio of the core 
and the neutron (4) and as well as on the nature of the 
subsystems, bound or virtual 0 The functional form 
of each observable written in terms of dimensionless ra¬ 
tios is given by a universal scaling function valid for short 
ranged interactions or large systems [l^. The model is 
more general and can describe phenomena in atomic and 
molecular physics such that other large three-body sys¬ 
tems with different structures can also be treated and 
stability and size studies performed (see e.g. Ref. 

We show that the asymptotic limit ^ I is reached 
in our model at much higher values of than the ones 
found in previous data analysis 0, B S 01 ■ Due to co¬ 
herence of the neutrons in the halo and final state inter¬ 
action, Cnn goes smoothly to the asymptotic limit only 
after displaying a minimum, as it will be shown. 

Particle distributions in the halo were calculated in 
Refs, mill, where the three-body Faddeev equations 
for the renormalized zero-range two-body interactions 
were solved. We obtain C„„ by using the correspond¬ 
ing three-body wave function with the inclusion of the 
nn final-state interaction (FSI) [^. We will show our 
results in case of the halo nuclei ®He, ^^Li and ^"‘Be. 

For the n — n — A three-body system, Cnn is given by 


Cnn — C {jpnn^ — 


nn l^(g nniPnn) \ 

fd^q 

nn P{qnA)p{qn'A) ’ 


( 1 ) 


j _^ ^nn 1 * —^ ^nn 

QnA — Pnn and Qn'A — Pnn 2 ’ 

where the one-body density is 


P{qnA) = / d^Qn'A 


-quA - qn'A, 


QnA Qn'A 


•( 2 ) 


$ = ^(qnn,Pnn) is the Corresponding breakup amplitude 
of three-body wave function including the FSI between 
the neutrons, qnn is the relative momentum between the 
core A and the center-of-mass of the nn subsystem; and 
Pnn the relative momentum between the neutrons. 

The FSI is introduced directly in the inner product 
$ = {qnn;Pnn^~'’\'^), where the ket |pn„ refers to 
the nn scattered wave given by the Lippmann-Schwinger 
equation. The correlation function, calculated with the 
distorted-wave amplitude, assumes a sudden breakup of 
the halo as the main reaction mechanism. The halo is 
considered as a coherent source of neutrons, differently 
from the framework of Ref. 0 In our picture, the slow 
halo motion decouples from the fast motion of the core 
in the field of the target. The distorted wave amplitude 
$ is given by 


^ — ^iqnnjPnn') T 


I/(27r2) 


V^rm '^Pn n 


Pin -p^ 


where is the three-body wave function 0 Cnn is the 
nn virtual state energy taken as 0.143 MeV. 

In the framework of Lednicky-Lyuboshits 0 , in order 
to obtain the correlation function Cnn, the probability 
density (in configuration space) of the neutron-neutron 
scattering state is multiplied by the probability density of 
the relative motion of the halo neutrons in the three-body 
wave function. Such model for the correlation function 
was developed from astrophysics, and the possibility to 
apply it to particle and nuclear physics is valid if one 
assumes that the particles emitted by the source are in¬ 
dependent. We believe this is not the case, when con¬ 
sidering the neutrons of the halo of a nuclei. So, in our 
approach, when the distorted-wave amplitude of Eq. © 
is considered, |<I>p contains off-diagonal matrix elements 
of the two-body densities. 

The results of our calculations for Cnn of the systems 
^^Be, ^^Li and ®He are respectively shown in Figs.^toOl 
They are shown as functions ofpnn, and compared with 
experimental available data B B 0 ■ 

In the present Faddeev renormalized zero-range model, 
as expected, Cnn goes asymptotically to unity, as a func¬ 
tion of Pnn- Such limit is shown clearly in all the three 
selected systems that we have analyzed. The results that 
include the asymptotic limit are shown in the insets of 
the figures. As we observe, for example, in the insets of 
Figs, n to El the theoretical asymptotic limits for 
the systems ^^Be, ^^Li, and ®He are, respectively, > 400 
MeV/c, > 250 MeV/c, and > 500 MeV/c. Besides the 
fact of the qualitative similar behavior, when comparing 
the results of the insets with the corresponding experi¬ 
mental results, we note a clear scale discrepancy. The 
striking observation is that the interference effect pro¬ 
duced by the inclusion of FSI originates a minimum for 
Cnn ~ 0.35, pushing the asymptotic limit to much larger 
values of Pnn than the ones considered in the asymptotic 
normalization of the experimental data. Our model calls 
for a different normalization of the experimental C'nn ob¬ 
tained in Refs. BB0- We believe that the qualitative 
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FIG. 1: Two-neutron correlation for the halo of ^'*Be, as a 
function of the relative nn momentum, Pnn- The solid curve 
gives the model results for S 2 n = 1.337 MeV, and EnA=0.2 
MeV. When compared with data, the model result is multi¬ 
plied by 1.425. Experimental data are from (open trian¬ 
gles) and (full circles). 


picture presented by our model, evidencing a minimum of 
Cnm will survive in a more realistic three-body approach. 

In Figs, n to 13 we observe a quite good agreement 
between theory and data when an appropriate normal¬ 
ization factor is included. Particularly, for ^^Be, the 
data are well reproduced by the model with S' 2 „=1.337 
MeV 1^ and with EnA =0.2 MeV (where EnAis the 
virtual state energy of n—^^Be), when our C„n is multi¬ 
plied by 1.425, as shown in Fig. ^ (or, Cnn ~ 0.70 
where is the model fit to experimental data). 

We have studied the model sensitivity to the two and 
three-body binding energies considering the case of ^^Li. 
We chose this nucleus because it has the smallest 2n sep¬ 
aration energy 821 % among the nuclei we are analyzing. 
The model results are shown by the three plots presented 
in the inset of Fig. [3 For a fixed s—wave ^°Li two-body 
virtual state energy, EnA =0.05 MeV 0, we vary S 2 n 
from 0.29 MeV (solid line) to 0.37 MeV [0 (dot¬ 
ted line). They differ only near the origin. Next, to ob¬ 
serve the sensitivity to the s—wave ^°Li two-body virtual 
state energy, we also plot a curve for EnA =0.8 MeV 0 
and S 2 n = 0.37 MeV (dashed line). Increasing the vir¬ 
tual state energy, the “plateau” near the minimum is 
enlarged and, for < 60 MeV/c, Cnn is enhanced and 
even strongly near zero. 

We should note that, as we increase the virtual state 
energy (reduce the absolute value of n—®Li scattering 
length) for a fixed S 2 n in ^^Li, we are shrinking the three- 


FIG. 2: Two-neutron correlation for the halo of ^^Li, as 
a function of the relative nn momentum, Pnn- The model 
results (inset) are given for three cases: S 2 n = 0.29 MeV and 
E„a= 0.05 MeV (solid line); S 2 n = 0.37 MeV and £'„a= 0.8 
MeV (dashed line); and, S 2 n = 0.37 MeV and EnA=0.05 MeV 
(dotted line). In the main body of the hgure, the solid curve 
(when r™® = 8.5fm fl3 i presents the corresponding curve 
of the inset multiplied by 2.5; the dot-dot-dashed curve, the 
model presented in 0 with r™® = 8.3 fm. The experimental 
data are from Q (full circles) and (empty circles). 


body system (see discussion and Table 1 of Ref. 0 ), 
implying that Cnn develops a longer tail in momentum 
space. In the inset of Fig. 13 this effect is clear and shown 
by the variation of the n—®Li virtual state energy from 
0.05 to 0.8 MeV, while S 2 n is kept fixed at 0.37 MeV. The 
increase of S 2 n from 0.29 to 0.37 MeV, with i5„Afixed at 
0.05 MeV, also implies in reducing the size of the three- 
body system. This is shown in the inset of Fig. 13 by com¬ 
paring the solid with dotted lines; the variation is small 
but consistent with the small variation of S 2 n- We ex¬ 
pect that such typical behavior is general. With respect 
to the correct asymptotic behavior, one can observe our 
solid line given in the inset figure. We need to multiply 
it by a factor 2.5 in order to fit the experimental data, as 
shown in the main body of the figure. So, in this case we 
have Cnn ~ 0.40 C!^. We also compare our calculation 
with another model, given in Q, which is close to the 
one presented in in the framework of 0 . We observe 
that the different behavior between our and other mod¬ 
els, essential for the absolute normalization of the data, 
will appear at higher momentum values. 

Next, to study the case of ®He, it’s worthwhile to men¬ 
tion that the s—wave n—'*He scattering has a positive 
scattering length and effective range , which by 
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FIG. 3: Two-neutron correlation for the halo of ®He, as a 
function of the relative nn momentum, Pnn- The solid curve 
gives the model results for S 2 n = 0.973 MeV and EnA=4: MeV; 
and, the dashed curve, for S 2 n = 0.973 MeV and E„a=0. 
When compared with data (extracted from |^), the model 
results (in the main body) are conveniently normalized as 
explained in the text. 


extrapolation of the standard effective range expansion 
produces an unphysical ®He bound state. Similar situa¬ 
tion is found in the quartet neutron-deuteron scattering 
where the large positive scattering length and an effective 
range produces, by naive extrapolation, a bound state 
near the threshold. The effective interaction in the case 
of the n — d system is repulsive and no quartet trinu¬ 
cleon exists, similarly for the n—^He where the s—wave 
interaction should be effectively repulsive to avoid an un¬ 
physical bound state. However, in our model the sign of 
the n—core scattering length is determinant of the occur¬ 
rence of the bound state. 

The experimental n—^He scattering length is 
3.26(3) fm 1^, which by naive extrapolation of 
the effective range in our zero-range model produces an 
unphysical bound state. Therefore, we arbitrarily made 
a calculation of Cnn with a zero energy s—wave ®He and 
the experimental value of S 2 n in ®He. Surprisingly, we 
found a reasonable agreement with the data, as shown 
in Fig. 121 In this case, our results for C„„ need to be 
multiplied by 1.48 to fit the data, as shown by the dashed 
line in the main part of the figure {Cnn ~ 0.68 C^). 


Observe that the interaction for the n—^He should be 
weaker than the one we are using; and, as the interac¬ 
tion becomes weaker while the three-body energy is kept 
constant, the three-body system tends to be more com¬ 
pact M implying in an increase of the typical three- 
body momentum scale. To mimic this effect, we decrease 
the attraction of the n—^He s—wave interaction, allowing 
a virtual s—wave state (unphysical), using as a typical 
magnitude for the scattering length the value of 3 fm. 
Still in this case, as we have multiplied our by 1.12, 
we obtain a good fit to data, as shown by the solid line 
of Fig.3. So, leaving the normalization free, one cannot 
distinguish between different weakly attractive s—wave 
interactions in the ®He, if the S 2 n and the nn virtual 
state energy are fixed. 

Alternatively, one may think in using the zero-range 
T—matrix and project out the bound-state, using the 
experimental value of the n—^He scattering length. This 
may effectively weaken the interaction, which could in 
principle produce a similar effect as we have discussed. 

In conclusion, we presented a theoretical study of 
neutron-neutron correlation for Borromean three-body 
weakly-bound systems, considering three cases where ex¬ 
perimental data are available. As the present analysis of 
halo dissociation in a renormalized zero-range three-body 
framework contains the main physical scales of weakly 
bound Borromean systems, we believe that our results 
for low nn relative energies (below 40 MeV) are to a large 
extend universal; i.e., other short-ranged potential mod¬ 
els, with the same low energy scales, will produce similar 
results. Considering that in our model there are no free 
parameters, as the inputs are just the physical scales (two 
and three body observables), and Cnn{Pnn) !> when 
Pnn oo, one can use it to get the normalization of the 
data. Our zero-range calculations suggest that the rele¬ 
vant effect of the inclusion of final-state interactions with 
the halo considered as a coherent source of neutrons is 
the presence of a minimum in Cnn at an intermediate 
region of (see Figs. lll3ll . It will be important to ver¬ 
ify this conclusion, as well as possible deviations, within 
a three-body model with more realistic two-body inter¬ 
actions. Finally, we observe a good agreement between 
theory and experimental results, provided that the data 
fits, C^, are appropriately normalized. In view of our 
findings, it will be very interesting to improve the ex¬ 
periments iiilli/ particularly at higher momenta, 
in order to characterize the existence of the minimum in 
Cnn- 
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